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Exotic axions 
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We show that axion phenomenology may be significantly different than conventionally assumed 
in theories which exhibit late phase transitions (below the QCD scale). In such theories one can find 
multiple pseudoscalars with axion-like couplings to matter, including a string scale axion, whose 
decay constant far exceeds the conventional cosmological bound. Such theories have several dark 
matter candidates. 
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I. INTRODUCTION 

The smallness of the experimentally determined upper 
bound on the strong CP violating parameter, 9 < 10~ 9 , 
is an outstanding puzzle of the standard model. One 
can either assume CP to be an exact symmetry spon- 
taneously broken in such a way as to ensure that 6 is 
naturally small, as in the Nelson-Barr mechanism 0, , 
or one can introduce a U(X) symmetry, known as the 
Peccei-Quinn (PQ) symmetry [3,3 to allow 8 to dynam- 
ically relax to zero. (Models suppressing strong CP vi- 
olation via spontaneously broken parity symmetry have 
also been constructed [H,|(|-) An attractive feature of the 
PQ mechanism is that it divorces the strong CP prob- 
lem from flavor physics — the masses and mixings of the 
quarks — whose origin remains a mystery. The PQ mech- 
anism entails a global U(l) symmetry which is exact 
up to a QCD (and possibly electromagnetic) anomaly. 
The symmetry breaks spontaneously at a scale /, giv- 
ing rise to a pseudoscalar Goldstone boson, the axion 
which couples to matter via the interaction (a/f)GG 
HH H, E3 El 13. Here a is the axion, / is its decay 
constant, G^ is the gluon field strength, and the ratio 
(a/ f) should be thought of as an angle. This angle has a 
potential arising from instantons which causes (a//) to 
sit at the vacuum where 9 — 0. The axion will in gen- 
eral have additional derivative couplings to matter, such 
as a model dependent coupling to photons of the form 
(a/f)FF [lj, \XM- The mass of the axion m a satisfies 
m a w m-xjir/f where m w ss 140 MeV and f n w 93 MeV 
are the pion mass and decay constant respectively. 

The axion decay constant is bounded from below by 
collider experiments and astrophysical arguments. The 
latter are the more stringent: if / is too small, the cou- 
pling to ordinary matter is large enough to allow copious 
axion production in red giants and supernovae, leading to 
an unacceptably rapid cooling rate. This yields the lower 
bound 10 9 GeV < / EH- It has also been shown that in 
conventional axion models, large / leads to copious pro- 
duction of cold, degenerate axions in the early universe 
El EEi, so that / > 10 12 GeV leads to unacceptably 
large f2d m , where fldm is the fraction of dark matter in 
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the universe today, measured to be 0.21 to within 4% [H^ |. 
Therefore the axion decay constant is conventionally as- 
sumed to lie in the window 10 9 GeV < / < 10 12 GeV. 
At the upper end of this bound axions are a viable dark 
matter candidate, and experimental attempts to detect 
their presence are in progress |2(j ■ 

In this Letter we will present a scenario for axions 
which exhibits exotic axion cosmology, and so we sum- 
marize here the conventional picture. The PQ symmetry 
breaks spontaneously at a temperature T ~ / well above 
the QCD scale where the instanton induced axion po- 
tential turns on. Following this phase transition (a//) 
equals some random angle 6i until T ~ 1 GeV, when 
the axion potential develops rapidly. The axion field be- 
gins to oscillate at temperature T$ ~ 1 GeV when the 
axion mass comes within the horizon, m (T,-) f=a H(Ti), 
where H(T{) is the Hubble parameter at temperature 
Tj, which is only weakly dependent on /. The coher- 
ent oscillation may be thought of as a gas of degenerate 
nonrelativistic axions with number density of axions per 
comoving volume equal to n a = 6j L H(Ti)f 2 . This quan- 
tity remains constant since the axions are nonrelativistic 
and because their annihilation rates are negligible. As a 
result, the subsequent axion energy density at tempera- 
ture T is given by p a = m a n a (R(Ti)/R(T) f, where R(T) 
is the Robertson- Walker scale factor at temperature T. 
The upper bound on / follows from observational limits 
on £ldm since p a varies almost linearly with /, assuming 
that 9i = 0(1). 

There have been prior attempts to evade the cosmolog- 
ical bound, motivated in part by the fact that in string 
theories, axions with / = y/2auM p w 10 16 GeV, where 
olxj is the unified value of the fine structure constant, are 
generic. A trivial way to harmlessly incorporate an axion 
A with decay constant F > 10 12 GeV is to introduce a 
second pseudoscalar a so that the Lagrangian contains 
the term (A/F + a/ f)GG. Provided that the simultane- 
ous transformation A —* (A+eF), a — > (a— ef) is an exact 
symmetry of the full theory, the spectrum will consist of a 
innocuous, massless Goldstone boson plus an axion with 
decay constant / subject to the usual constraints sum- 
marized above. Another resolution is to assume inflation 
and an ensemble of initial angles 9i to choose from, in- 
voking the anthropic principle to justify a value in our 
universe of Oi s=s |2^|; recently it was pointed out that 
such a scenario could be constrained by the PLANCK po- 
larimetry experiment |22| . In addition, Banks, Dine and 
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Graesser have pointed out that any attempt to evade the 
cosmological upper bound on the axion energy must deal 
with the generally more difficult question of the cosmo- 
logical saxion energy problem |23j |. Another interesting 
proposal by Barr and Kyae (BK) |24| is an axion cou- 
pled to the light quarks in such a way to allow for ex- 
otic cosmological evolution of the axion below the QCD 
scale. Our motivation here is to explore unconventional 
axion cosmology without invoking the anthropic princi- 
ple, and without convoluting the PQ mechanism with 
flavor physics; in common with the BK suggestion, and 
the earlier ref. |25| . our scenario involves late evolution 
of the PQ symmetry breaking order parameter, although 
it differs significantly in realization and phenomenology. 

The conventional cosmological bound on / results be- 
cause the PQ symmetry spontaneously breaks at a tem- 
perature well before the QCD scale. However, if / only 
evolves to a large vev after the QCD time, such bounds 
may be evaded. This could be accomplished with a suf- 
ficiently flat PQ potential, so that the radial mode only 
evolves out to a large vev once its mass enters the hori- 
zon. 

The fine tuning associated with such a flat potential 
can only be avoided with supersymmetry (SUSY); but 
even then SUSY breaking will in general generate a cur- 
vature for the PQ potential which forces PQ breaking at 
or above the weak scale, a transition which is too early 
to evade the cosmological constraints. This problem can 
be avoided with the introduction of a new sector which 
couples only weakly to the standard model via the con- 
ventional PQ sector. We show that such a coupling may 
be weak enough to shield the potential in the new sec- 
tor from SUSY breaking effects, while still significantly 
affecting axion phenomenology. The spectrum of this 
theory includes (i) an axion far lighter than convention- 
ally allowed with no significant cosmological abundance; 
(ii) an additional axion-like particle which is heavier than 
would be an axion with comparable decay constant; (iii) a 
light dilaton-like scalar particle. The latter two are dark 
matter candidates. In the next sections, we describe a 
model which realizes this scenario. 



II. A MODEL 

Our starting point is to assume a viable supersymmet- 
ric theory which implements the conventional PQ mech- 
anism. We assume that there exists a superfield <j>i which 
carries PQ charge and couples to colored fermions in a 
real representation of the gauge group; at a tempera- 
ture well above the QCD scale this field acquires a vev 
(4>i) = vi/y/2 which lies within the conventional window 
10 9 GeV < vi < 10 12 GeV. With this vev, the heavy col- 
ored fermions coupled to (f>i develop a mass Mq — gv\. 
It is important for our modification of the theory that 
that the saxion be light (to be specified below) so that it 
not communicate large SUSY breaking to a new sector 
we will be adding. A light saxion is expected in any the- 



ory of low energy SUSY breaking; it could also occur in 
gravity mediated SUSY breaking models so long as the 
PQ sector is sequestered from the SUSY breaking. An 
excess of saxion energy can be avoided in such models 
either by having relatively late inflation (below the PQ 
breaking scale but well above the QCD time), or by hav- 
ing the minimum of the saxion potential be at the same 
point as preferred by finite temperature effects prior to 
an epoch of higher scale inflation. 

To this theory we now introduce the supcrpotential 

W = V2XA(h^4> 2 -4>l) + ^ 2 B (2(<VV) 2 - 1) » (1) 

where 4>2 carries opposite PQ charge from (f>i, while A,B 
and (j>o are PQ-invariant fields. The four parameters Vq, 
fi 2 , h and A may all be taken to be real by redefinition 
of the phases of A, B, (f> and (j>2- The part of the scalar 
potential relevant to us is 

Vfa, <h, <M = 2A 2 |/^ 2 - + ^\2(cb /v ) 2 - 1| 2 (2) 

The first term in V exhibits a flat direction in (f>2 and (f>o, 
which is lifted slightly by the second term in V, with A <C 
1 and /i -C Aqcb. The smallness of the couplings leave 
V / 2(0i) = v\ unaffected, and the minimum of the almost 
fiat direction is at y/2((j>o) = Vq ^> V\, and y / 2(02) = 

V2 = vl/(hvi) > Vq. 

We assume that following inflation and reheating the 
Universe sits away from the minimum of the potential 
with V2(4>i) = «i, and (A) = (B) = (fa) = (<fe) = 0, 
the latter being determined by high temperature effects 
due to interactions with unspecified heavy fields prior 
to inflation. This field configuration persists down to 
a temperature To < Aqcd, satisfying ^ 2 /vq s» H(Tq), 
when the curvature lifting the flat direction is sufficiently 
strong to overcome the Hubble friction. Then </>o will 
roll out to its minimum at (<po) — vq 3> Wi, causing cf>2 
to follow its flat direction out to V2 = v 2 /(hvi) 3> vq. 
With [i <C Aqc d this phase transition will occur at Tq < 
Aqcd, and as we discuss further below, the dark matter 
produced in this transition can be made acceptable. 

To understand how the pseudoscalars behave in this 
model, we expand around the minimum of the potential, 
writing ^ = -^(vi + ai)e l7Ti / Vl with i = 0, 1, 2. The fields 
0% and TTi are the scalar and pseudoscalar excitations re- 
spectively. After adding to V the QCD contribution to 
the 7Ti potential, the complete pseudoscalar potential is 

T7 2,2 ""1 X 2 4 f 71 ! , ^2 27T \ 

= -m w f w cos A v cos 1 

Vl \«1 V 2 Vq J 



(up to an overall 0(1) factor 01 m front of the first 
term which does not concern us here). Since 7Ti is the 
only field that couples to ordinary matter, to understand 
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phenomenology we must decompose tt\ into mass eigen- 
states QSi,2,3- We diagonalize the mass matrix to leading 
nonzero order in 8 = vx/vq and e = p A / (X 2 v 4 ), and to 
all orders in x — A 2 Uq/(to^./^), finding the masses and 
decay constants 

m 2 ai = ^f-(l + x) , f 1=Vl , 
2 mlfl ( Ax \ f 1 + x \ ,a\ 

2 

TO a 3 = -2 . 
"1 

where the decay constants are defined by 

£qcd = —GG = f^y- + GG . (5) 
Vi V Ji J2 / 

In these formulas we have neglected terms of order vo jvi ; 
the 03 pseudoscalar decouples from the standard model 
to the order we work. 

This model yields a distinct particle spectrum: a\ is 
mostly 7Ti and couples like an axion with decay constant 
fx = Vx, except that its mass can be much heavier than 
that of a conventional axion by the factor y/l + x; ai 
looks like a very light axion, primarily composed of ttq, 
with decay constant f^ = Vo(l + x)/2x w vq/2. Since 
Vq 3> Vx, rn,2 is much lighter than that allowed by the 
conventional axion bounds — it is potentially an axion 
with string scale PQ constant f^ ~ 10 16 GeV. The en- 
ergy originally in nx is primarily transferred into the a\ 
field, making it a good dark matter candidate; the light 
a.2 axion receives only a small fraction of that energy, 
suppressed by (vx/vq) 2 -C 1. 

III. CONSTRAINTS 

At temperatures below the QCD scale but above the 
secondary transition, the universe has a background den- 
sity of cold 7Ti pseudoscalars which behave like a conven- 
tional axion with decay constant v\. Subsequently the 
(j>o field rolls out to its minimum ((j) ) — vq at the tem- 
perature T where H(T ) ~ p 2 /v . At this point (i) the 
energy in ttx pseudoscalars gets redistributed among the 
Oi and a 2 mass eigenstates following eq. @; (ii) an en- 
ergy density /i 4 is released and is mostly transferred into 
txo (radial) oscillations of the 0o field. There are two cos- 
mological constraints on the late transition. The first is 
that the combined energy density in the exotic ax axion 
and do scalar not exceed the observed dark matter den- 
sity today. A second constraint is that the phase transi- 
tion occur after the QCD time but before today (or well 
before matter-radiation equality if either ax or uq con- 
tribute an appreciable fraction of the dark matter). In 
addition there is the astrophysical constraint that to pre- 
vent copious production of any of the light pseudoscalars 
in supernovae, we require fx, fi ^ 10 9 GeV. We now 




10 9 10 10 10 11 10 12 

fl GeV 



FIG. 1: The mass of the exotic axion a\ vs. its decay con- 
stant fi. Two solid parallel lines give m ai (fi) assuming ai 
comprises half the dark matter for initial misalignment angle 
Oi =0.1 (upper line) and 9i = 1 (lower line). Dashed line is 
the axion mass in standard PQ models, with the heavy dashed 
line giving the region where conventional axions could be the 
dark matter. 



discuss these constraints in detail, and map out the cor- 
responding parameter space. 

In order to ensure no excess of dark matter, we must 
limit the energy in nx by requiring fx < 10 12 GeV. We 
also require that the energy po in the form of oo oscilla- 
tions produced at the secondary transition not dominate 
the universe at the epoch of matter-radiation equality, 
T eq sb 1 eV: 

Po{T eq )^^(^X <Tt q , ff(T )«^. (6) 

A complication arises from the fact that energy in the 
early ttx oscillations is transferred primarily into ax oscil- 
lations, while ax becomes heavier than nx by a factor of 
VI + x as 4>o rolls from 4>q = to </>o = ^o- This increase 
of energy must come from the energy released during the 
secondary transition; it can be represented by a contri- 
bution to the potential for <f>Q due to the (^-dependent 
energy density of the ax axion condensate |2^, of the 
form 



Px(0o,T) = n Vl (T) m ni ^/l + x(4> ) , (7) 

where Tj ~ 1 GeV is the temperature when ttx begins 
to oscillate, n Vl {T) = (OfHiff) (T/Tj) 3 is the number 
density of nx bosons in the temperature range Aqcd > 
T > To, TO^ = {rnirfn/ fx) is the tti mass before <j>o rolls 
out, and x(<f>o) = (4>q/vq) A x controls how the tti mass 
becomes the ax mass as <fio increases. If px(vq,Tq) > p 4 , 
this potential delays the secondary transition, so that <f>o 
only gains its vev at some lower temperature, T = Tq 
satisfying px(vo,Tq) w Po(Tq) rs p, 4 . After this delayed 
transition both the energy density in 4>q and a\ oscilla- 
tions remain comparable, diluting cx T~ 3 . After a little 
algebra one finds that eq. © still holds in this case, 
but is augmented by the constraint that there not be too 
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TABLE I: Two sets of parameters (1,11) allowed by the con- 
straints eqs. H6I8I101 . Set I gives parameters for an prompt 
transition at temperature To, set II gives parameters for a 
delayed transition at Tq. All parameters are in GeV except 
the axion masses which are in eV. 

much energy in the cosmological a\ abundance, 

Tt q > (T^ n ^{T Q )m ai , (8) 

with m ai defined in eq. and To ~ W M p p? /vq from 
cq. If we now assume pi+ po = pdm, that is that ai 
and do together compose the dark matter, we can use eq. 
||SJ| to compute allowed masses and couplings, /i, m ai , 
for the observable dark matter candidate a±. The result 
is shown in fig. 1; evidently the mass and coupling of 
a\ can differ significantly from those of the conventional 
axion. 

Finally, we obtain a constraint from requiring fa to 
roll after the QCD phase transition, but somewhat before 
matter-radiation equilibrium if either a\ or <jq are to be 
the dark matter: 

/ \ 1/2 / \ 1/2 

If the origin of dark matter lies elsewhere, then T eq in 
the above equation is replaced by today's temperature. 

The model as it stands possesses an exact discrete sym- 
metry S — > —S which is spontaneously broken and leads 
to domain walls. This can be avoided by breaking the 
symmetry explicitly, either with a small linear term in 
the low energy superpotential eq. (JIJ, or in the high 
energy interactions such that (S) ^ (but < v^) after 
post-inflationary reheating. 

Aside from cosmological constraints, there is also a nat- 
uralness constraint. Having a flat potential is critical for 



the late phase transition, so we require that SUSY break- 
ing terms induced by the interactions between fa , fa and 
fa do not lift the flat direction too much. The most im- 
portant term is the soft mass generated for fa from the 
SUSY breaking fa mass, for which we require 

2 2 

A 2,2 m a 1 m < p 1 ^4 fm\ 
2/2 1&J ~ ■ (10) 

This imposes a significant new constraint on p. In ad- 
dition, we are assuming this sector is sequestered from 
gravity- mediated SUSY breaking pp| . 

The combination of constraints eqs. (|6I8I10|I yields a 
parameter space too large to explore in detail here [2^ |. 
We give here instead two representative sets of values in 
Table I satisfying the constraints. 

IV. PHENOMENOLOGY 

Conventional axion models have fairly circumscribed 
phenomenology; if one assumes that the dark matter con- 
sists of a conventional axion, then the mass and coupling 
of the axion are related in a direct way, and both lie in 
a fairly model independent range about m a ~ 10 -5 eV, 
f a ~ 10 12 GeV determined by the initial axion misalign- 
ment 6i. If the axion is not the dark matter, its mass can 
be heavier. 

In contrast we have shown how these simple relations 
can be greatly modified in a theory with a late phase 
transition below the QCD scale, as shown in table I. Dark 
matter in this theory can be comprised of roughly equal 
parts of remnant pseudoscalar and scalar particles. 

It will be interesting to examine whether the late phase 
transition discussed in this Letter could give rise to iden- 
tifiable signatures, due to large isocurvature fluctuations 
seeded at the time of the transition, over angular scales 
which are small relative to the measured microwave back- 
ground fluctuations. 

We are grateful for numerous conversations about this 
work with C. Hogan and A. Nelson. This work was sup- 
ported by DOE grant DE-FGO3-00ER41132. 
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